Harmonic forcing of an extended oscillatory system: 
Homogeneous and periodic solutions 
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In this paper we study the effect of external harmonic forcing on a one- dimensional oscillatory 
system described by the complex Ginzburg-Landau equation (CGLE). For a sufficiently large forcing 
amplitude, a homogeneous state with no spatial structure is observed. The state becomes unstable to 
a spatially periodic "stripe" state via a supercritical bifurcation as the forcing amplitude decreases. 
An approximate phase equation is derived, and an analytic solution for the stripe state is obtained, 
through which the asymmetric behavior of the stability border of the state is explained. The phase 
equation, in particular the analytic solution, is found to be very useful in understanding the stability 
borders of the homogeneous and stripe states of the forced CGLE. 

PACS numbers; 05.45.Xt, 89.75.Kd, 47.20.Ky 



I. INTRODUCTION 



Nonequilibrium pattern formation is widely observed 
in many physical, chemical and biological systems. Sig- 
nificant progresses have been made in this field during 
the last few decades. For example, it has been found 
that nonequilibrium patterns can be grouped into a few 
universality classes In many cases, such a system 

is in constant interaction with its environment, and un- 
derstanding the effect of extrinsic perturbation is of great 
theoretical and practical importance. In particular, it is 
interesting to study the deformation of an existing pat- 
tern or the formation of a new pattern under an external 
forcing. Our understanding in such a direction is far from 
complete. 

Petrov et ai, and later Lin et ai, studied the light 
sensitive Belousov-Zhabotinsky (BZ) reaction in an os- 
cillatory regime in the presence of a periodic modula- 
tion of the intensity of illumination [||J|]. They ob- 
served "entrainment bands" in which the system is fre- 
quency locked. Different spatial patterns — stationary 
fronts, standing waves of labyrinth and more complex 
shapes — are observed within the bands. In a similar BZ 
reaction setup, Vanag et at studied spatial patterns and 
transitions among them in detail, and observed localized 
irregular/standing clusters as well as the above patterns 



Continuum models of forced pattern forming systems 
can be grouped into ones based on a kinetic model or on 
an amplitude equation. In the first group, an unforced 
system is modeled by a coupled kinetic model, such as the 
Brusselator, Oregonator or FitzHugh-Nagumo model, 
and parameters in the model are modulated to simulate 
the effect of external forcing (e.g. |^,^). On the other 



hand, near a bifurcation onset of a pattern, small differ- 
ences among systems become irrelevant, and they are all 
described by one of a few universal equations. If the bi- 
furcation is supercritical and oscillatory, and if the most 
unstable wave number is zero, the complex Ginzburg- 
Landau equation (CGLE) is the equation for the class of 
systems. In the presence of an external periodic modula- 
tion, it is shown that the CGLE with an additional forc- 
ing term becomes the appropriate equation |^,^. There 
exist a few studies on forced CGLE, and diverse behav- 
iors are observed depending on several factors, such as 
the spatial dimension, the mode of the frequency locking, 
and the behavior of the corresponding unforced system 
^,P|-[l6| . However, we do not even know what behaviors 
are possible, let alone understand them. 

Even the simplest case of the 1:1 locking in one dimen- 
sion displays a wide variety of behaviors. At large ampli- 
tude of the forcing, a homogeneous state with no spatial 
structure is stable. Chate et al. found that the homo- 
geneous state becomes unstable to a periodic "stripe" or 
"kink-breeding" state as the forcing amplitude decreases, 
and a "turbulent synchronized" state — chaotic with its 
average phase is locked to that of the forcing — can ap- 
pear, as the amplitude decreases further fl^ . 

In this paper, we study in detail the homogeneous and 
stripe states of one dimensional forced CGLE around the 
1:1 locking. There are two borders regarding the homo- 
geneous state (1) the stability border, below which the 
state loses its stability, and (2) the existence border, be- 
low which a homogeneous solution does not exist. In gen- 
eral, the existence and stability borders do not coincide. 
It is known that the stability border of the homogeneous 
state lacks a reflection symmetry around the v — a line. 
Here, v is the difference between the natural and exter- 
nal frequencies, and a is a nonlinear detuning parameter. 



1 



We find the asymmetry can be explained by the linear 
stability of the state. Also, the condition under which 
the existence and stability borders of the homogeneous 
state coincide is found. The stability border of the stripe 
state also lacks a reflection symmetry. An approximate 
phase equation is derived from the forced CGLE, and it 
is found that its qualitative behavior is identical to that 
of the original equation, at least in the region of present 
interest. An analytic expression of the stripe state for 
the phase equation is obtained, which is used to explain 
the asymmetry of the border of the stripe state. 



II. FORCED COMPLEX GINZBURG-LANDAU 
EQUATION 

A. Complex Ginzburg-Landau equation 

Near the stability border of a homogeneous state of an 
extended pattern forming system, the time evolution in 
a large spatial and temporal scale is given by one of a 
few universal equations . If the instability is oscilla- 
tory and supercritical, and the wave number of the most 
unstable mode is zero, the CGLE, 



dtA = A-{\ + ia)\A\'^A + (1 + iP)V'^A 



(1) 



is the governing equation. Here, A is complex amplitude, 
and a, (3 are real constants. The behavior of the CGLE is 
relatively well understood, especially in one and two di- 
mensions |r^-|l9|]. It has plane wave solutions, which are 
stable only if 1 -I- a[3 > 0. Otherwise, the Benjamin-Feir 
instability sets in, making the solutions unstable. Near 
the unstable side of the stability border {l + a/3 — line), 
"phase turbulence" is observed, which is characterized by 
disordered cellular structure and the absence of a defect 
(1^1 = 0). "Defect turbulence" is observed further in 
the unstable region, where constant creation and annihi- 
lation of defects is observed [^^. In this paper, the 
value of a = —3/4, /3 = 2 is mainly used, which is in the 
phase turbulence region. 



B. Homogeneous state 

Consider the case wherein a sinusoidal forcing is ap- 
plied to the system of Eq. (|l|). It was shown that an 
additional forcing term should be included, and its form 
can be determined from the conditions of the spatial and 
temporal translation invariance . For a harmonic forc- 
ing (near the 1:1 tongue), the resulting equation is 

dtA = (l-f ii^)A - (1 + ia) \Af A+{1 + iP)V^A + B, 



where u is the difference between the natural and forc- 
ing frequencies, and B is related to the amplitude of the 
forcing. 

We first seek for the homogeneous solution of Eq. (||) . 
In polar coordinates [A = Rexp{i^)], the equation be- 
comes 

dtR = R- + Bcos^ 

+R^, - /3i?$,, - 2/3i?,$, - R<^>1, 
Rdt^ = vR- aR^ - B sin $ 

+(3Rx. + + 2i?, - f3R^l . (3) 

For a sufficiently large B, the system is expected to lock 
to the forcing. Its homogeneous solution is 



Bcos$o = -i?o(l - Rl): 
Bsin$o = Ro{i^ ~ oiR^), 



(4) 



which can have one or three roots depending on the pa- 
rameters. For the three root case, only the one corre- 
sponding to the largest i?o is stable. The region of the 
parameter space in which a locked homogeneous solution 
exists is shown in Fig. 0(a). 

We apply the linear stability analysis to the homoge- 
neous solution ||l^, where the behavior of a small devi- 
ation from the solution r = R — Rq and = $ — $o is 
studied. The growth rate of the mode with wave number 
k is found to be 

X{k) ^l-2Rl-k^ 

+ + a^)Rt -{v- 2aR^ - /3A:2)2, (5) 

which has the maximum value of 

A„,ax = 1 - 2i?2 _ - {2a + v/(l + a2)(l + /32)}i?2] 

(6) 

at fc = fcinax, corresponding to the most unstable mode, 
which is given by 



2 

max 




(7) 



(2) 



The stability border Bs of the homogeneous solution is 
obtained by solving numerically Amax — 0, which is also 
shown in Fig. |l|(a). 

A distinct feature of the stability border is that it is not 
symmetric to the v — a line. As shown in the figure, the 
difference between the existence and the stability border 
is smaller at the ly < a side. Moreover, the difference van- 
ishes for v < Vc with i^c — —1.067. This feature can be 
understood from the v dependence of femaxj which is given 
by Eq. (|7|) . It is found that femax is an increasing function 
of I/: it is zero for <i'c = [2a+ ^(1 + q;2)/(1 + /3)]i?g, 
and proportional to — slightly above v^. Since 
the wave number of the most unstable mode is zero for 
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v < Vc, and since the solution Eq. (Q) with the largest Rq 
is stable to a zero wave number perturbation, the exis- 
tence of the homogeneous solution guarantees its stabil- 
ity. 

As will be discussed later, an approximate phase equa- 
tion is derived from Eq. (||), which gives an additional 
insight into the stability border. The origin of the in- 
stability of the homogeneous state of the phase equation 
can be traced to a Laplacian term, whose coefficient is 
a decreasing function of i/, and becomes negative at i'^. 
Thus the homogeneous state is stable for v < vf , and it 
becomes more unstable as v increases. 



C. Stripe state 

The behavior below the stability border is investigated 
numerically. The forced CGLE in one dimension is inte- 
grated using a pseudospectral method for various v and 
B 1^ . The spatial resolution Aa; and time step Ai used 
are 0.1 and 0.01, respectively. Also, a periodic bound- 
ary condition is used. For most cases, the linear size of 
the system is chosen to be 4096, and the time interval of 
2 X lO** is used. Larger systems for longer intervals are 
also studied, and no change in the behavior is observed. 

The numerical integrations confirm the prediction that 
the homogeneous state is stable above the stability bor- 
der. It is found that the state undergoes a supercritical 
bifurcation to a spatially periodic static "stripe" state as 
B decreases below the border, and the modulation am- 
plitude of the state, defined as 54) = {{(f) — {(t>)xY)x, 
behaves as y/Bg — B close to the border [Fig. ||(a)] . The 
wave number of the stripe state is found to agree very 
well with fcmax of Eq. (0), especially near the border. 
In order to check how the nature of the transition de- 
pends on the unforced dynamics, the transition from a 
homogeneous state is examined for four different values 
of (a,/3): (-2,2), (-1.11,1), (-2,0), and (-0.75,0.5). 
It is found that a supercritical transition to the stripe 
state is observed for the first two cases belonging to the 
Benjamin-Feir (BF) unstable region, while a transition 
to a disordered structure is observed for the other two 
cases belonging to the BF stable region. 

As B decreases further, the stripe state becomes un- 
stable to a fluctuating stripe or "kink-breeding" state, 
depending on v p^ ]. The stability border of the stripe 
state is determined, and is plotted in Fig. |l|(a). Again, 
the border is not symmetric to the u = a line. The re- 
gion of the stripe state is much broader on the v > a side. 
Moreover, it extends to the region where a locked homo- 
geneous state does not exist. The origin of the asym- 
metry will be discussed using a phase equation, and is 
found to be very different from the case of the homoge- 
neous state. 



III. PHASE EQUATION 

A. Derivation 

An approximate phase equation can be derived from 
Eq. (H) as follows. Define small variables r — R— Rq and 
(j> = $ — $01 and assume that the time scale for <j) is much 
larger than that for r. The variable r is then slaved to (j). 
Starting from Eq. (||), it can be shown that 



Rq 



(l-i?^) + — COS($0 + ' 
-n-O 



(8) 



where additional terms higher than the second order in (p 
are ignored, and B is assumed to be small. Substituting 
this to the phase part of Eq. (|^), 



Rodtcj) = vRq - aRl - B\/l + a'^ sin($o + (j) + 5) 

+b(j)l + ctpxx + d(f>xxxx + e, (9) 

where a,b,c,d,e,S are constants depending on a,(3,i), 
and Rq. Since i?o is 1 at = a, and is a slowly varying 
function of it is expected that setting Rq = 1 does not 
change the qualitative behavior of the equation. On the 
other hand, the constants are simplified to 

a = a — {i) — a) /2, 

6 = a- /3-(i/- a)/2, 
c= l + a/3-/3(i/-a)/2, 

e = 0, 

5 = tan"i(a). (10) 

For the remainder of the paper, Rq will be set to 1 in 
the equation. Note that Eq. (|^) is a generalized version 
of the phase equation obtained by CouUet and Emilsson, 
which is derived for the special case of ~ a |1. Also, 
4>xxxx term is added for the stability of the solution in 
the phase and defect turbulence regions. 



B. Homogeneous state 

The phase equation is studied in a way parallel to the 
analysis of the forced CGLE. Homogeneous states, given 
by 



$0 = sin ^{ 



bVT 



f)-S 



(11) 



exist for B > {v — a) / VT+a^ . There exist two solutions 



$0 and $g, in the [0, 27r] interval satisfying Eq. (11) as 
shown in Fig. |^. The $5 solution is stable under homoge- 
neous perturbation, while the $o solution is unstable. A 
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linear stability analysis of the stable homogeneous state 
shows that the maximum growth rate is 



-SVl + a^ cos($o + 5)- c^/M 



(12) 



for the mode with wave number k^^^ = y^c/2d (if c < 0). 
This state is found to be linearly stable above the stabil- 
ity border Bf, which is given as 



(cV4ci)^ + (i^ - ay 

rT^2 



The existence and stability borders are plotted in 
Fig. |l](b) for a = -3/4, /3 = 2. Note that the shapes 
of the borders are qualitatively the same as those of the 
forced CGLE: the stability border is asymmetric to the 
1/ — a line, and the two borders meet for v < vf . Since 
the wave number of the most unstable mode should be 



zero at v ■ 



one arrives at 



2 + 3a/3 



(14) 



For the above values of a and /3, vf — —5/4, which is 
comparable to the value of for the forced CGLE. 

The simple structure of the phase equation makes its 
interpretation simple. The reason for the instability is 
that c can be negative, while the (j^xxxx term always tries 
to suppress such an instability. The value of c remains 
positive for v < vf , and the homogeneous state is sta- 
ble. As V increases further, c becomes negative. Since 
c is a decreasing function of the instability becomes 
stronger with increasing which explains the fact that 
the difference between the existence and stability border 
increases with u. 



C. Stripe state 

The behavior of the phase equation below the stability 
border is studied numerically. As one crosses the bor- 
der, the homogeneous state goes through a supercritical 
bifurcation to a stripe state, and the modulation ampli- 
tude 5(j} = a/ {{4> — {(f>)xY)x behaves as yBt — B close 
to the border. A typical dependence of 50 on B is shown 
in Fig. ^(b), where = — 0.75. As B decreases further, 
the stripe state becomes unstable. The stability border 
of the stripe state determined numerically is plotted in 
Fig. 0(b). Again, the border is not symmetric to the 
V — a line, and even extends below the existence border 
of the homogeneous state. Although the phase equation 
is simpler than the forced CGLE, their qualitative behav- 
iors are essentially the same, at least for the homogeneous 
and stripe states. 

The simple structure of the phase equation allows an 
analytic expression for the stripe solution. The solution 
may be expanded in terms of harmonic functions 



(i}(x) = (pi+Si sin {kox) + Ci cos (kox) 

+S2 sin {2kox) + C2 cos {2kox), 



(15) 



where higher harmonics are ignored, fcp is the wave num- 
ber of the most unstable mode, and the coefficient Ci 
can always be set to by choosing an appropriate origin. 
Substituting it in Eq. (y), we find 



(13) sin($o- 



^2 

bi+S) = 



0, 



(5? 



-4C|) 



bc/M 



Sf - 4^2 



C. = 2 



bVTT^ (1 - {Sf + C|)/4) 
~BVl + q2 co8{<^o+(j)i + 5) + 2c^/d 

bVT 
-bVT 



o? sin(<i>o + 0i + (5) - bc/d 



bVT 



(16) 



a2 sin($o + </'i + '^) + 2bc/d 

which can be solved numerically for Si , C2 and 0i . Near 
the stability border Bf , an approximate analytic solution 
can be obtained, which is 



Sio^\jBt 
C2 oc Si, 



1 



, {v-a+^)Sl, (17) 

4B\/rT^cos($o + (5) ^ d' ^ ' 

where the proportionality constants are rather complex 
except for the case of 0i. The analytic solution agrees 
well with the results using numerical integration: as 
shown in Fig. ^(b), the modulation amplitude Scj) vs 
B curve obtained from the above expression is in good 
agreement with the corresponding numerical values. 

The analytic solution confirms not only the square 
root dependence of Ai, but it also provides an expla- 
nation for the asymmetry of the stability border of the 
stripe solution. Shown in Fig. ^(a) are two homo- 
geneous solutions — stable $q and unstable $g — of the 
phase equation with B a little below the stability border. 
As B decreases from the homogeneous toward the stripe 
region, the modulation amplitude around $0 increases 
with decreasing B. For sufficiently small B, (j){x) at cer- 
tain X approaches the unstable fixed point $q, which 
then makes the stripe solution unstable. Note that (pi 
in Eq. (^) is nonzero — it is negative when v is not very 
different from a. Thus, the average phase of a stripe 
state is shifted toward a value smaller than <i>Q. The av- 
erage phase of the stripe solution measured using numer- 
ical integration is also shown, which confirms the shift. 
As shown in Fig. ^, the shift moves the stripe solution 
toward (away from) <i>g for v < a > a), resulting 
in the asymmetry (a related argument is given in psf). 
The situation is entirely similar for the forced CGLE. 
The phase of the stable and unstable solution is plotted 
against v in Fig. ^(b) |2^. Also plotted is the average 
phase of the stripe solution. Here again, the average 
phase is shifted toward (away from) the unstable solu- 
tion for u < a > a). 
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IV. CONCLUSION 

Despite its simplicity, the forced CGLE displays a 
large variety of phenomena. The homogeneous and stripe 
states are mainly discussed here, and the phase equation 
is found to be very useful in understanding the stabil- 
ity borders of the forced CGLE. For a sufficiently large 
forcing amplitude, a homogeneous state with no spatial 
structure is observed. The state becomes unstable to a 
spatially periodic stripe state via a supercritical bifurca- 
tion as the forcing amplitude decreases. We obtained an 
analytic solution for the stripe state of the phase equa- 
tion, through which an argument for the asymmetry of 
the stability border of the state is formulated. As B de- 
creases further, more complex behaviors such as, the kink 
breeding, "depinning," and "roughening," are expected, 
which are currently under investigation. 
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FIG. 1. (a) The existence and stability borders of the ho- 
mogeneous state for the harmonically forced CGLE Eq. (^) 
with Q = —3/4 and /3 = 2. The stability border of the stripe 
state is also shown. Note that the stability borders are not 
symmetric to the v = ol line, (b) Corresponding borders for 
the phase equation Eq. (^ with a = —3/4, /3 = 2, and _Ro = 1. 
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FIG. 2. (a) Modulation amplitude 5(l> of the stripe solution 
of the forced CGLE Eq. is shown against B for f = —0.75. 
Also shown is a square root fit D\/ Ba — B with D = 1.25 and 
Ba = 0.02101. (b) 5(j> of the stripe solution of the phase equa- 
tion Eq. (^) is shown for v = —0.75. Analytic expression 
Eq. ( |l5| ) with the coefficients given by Eq. is found to be 
a good approximation. 
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FIG. 3. A schematic view of the determination of the ho- 
mogeneous solution of the phase equation Eq. (^: There exist 
one stable $o a^nd one unstable $o solutions. Here, the one 
with positive cos($o -f 5) is stable. For v > a,(^Q > while 
"I>o < ^0 for 1^ < a. 
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FIG. 4. (a) The stable $o sluA unstable $o homogeneous 
solutions of the phase equation Eq. (|^) just below the stabil- 
ity border {B — Bf = 5 x 10~^) are plotted against v. The 
average phase of the stripe solution of the equation for the 
same B is also plotted. The average phase is smaller than <l>o, 
which moves away from (toward) <l>o for v > a (v < a), (b) 
The same plot for the forced CGLE. The qualitative behavior 
is identical. 
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